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, Abstract. We have derived abundances of 33 elements and upper limits for 6 additional elements for the metal-poor ([Fe/H] = 

-2.42) turn-off star HE 0338-3945 from high-quality VLT-UVES spectra. The star is heavily enriched, by about a factor of 
100 relative to iron and the Sun, in the heavy i-elements (Ba, La, ..). It is also heavily enriched in Eu, which is generally 
considered an r-element, and in other similar elements. It is less enriched, by about a factor of 10, in the lighter .s-elements (Sr, 
Y and Zr). C is also strongly enhanced and, to a somewhat lesser degree, N and O. These abundance estimates are subject to 
severe uncertainties due to NLTE and thermal inhomogeneities which are not taken into detailed consideration. However, an 
interesting result, which is most probably robust in spite of these uncertainties, emerges: the abundances derived for this star 
i-Ch are very similar to those of other stars with an overall enhancement of all elements beyond the iron peak. 

We have defined criteria for this class of stars, r+s stars, and discuss nine different scenarios to explain their origin. None 
of these explanations is found to be entirely convincing. The most plausible hypotheses involve a binary system in which 
the primary component goes through its giant branch and asymptotic giant branch phases and produces CNO and i-elements 
which are dumped onto the observed star. Whether the r-element Eu is produced by supemovae before the star was formed 
(perhaps triggering the formation of a low-mass binary), by a companion as it explodes as a supernova (possibly triggered by 
mass transfer), or whether it is possibly produced in a high-neutron-density version of the i-process is still unclear. Several 
suggestions are made on how to clarify this situation. 
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1. Introduction 

Elements with atomic numbers Z > 30 are believed to be 
almost exclusively synthesized in neutron-capture (n-capture) 
reactions. In the most metal-poor stars the overall abundance 
of these elements varies from star to star, by more than a 
factor of 100 at a given metallicity (Mc William et al. 1995 
Ryan et al. 1996 1. Also, the different abundance ratios vary, 
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e.g. the Ba/Eu ratio tends to decline with decreasing [Fe/H] 
(Mc William 1998 Burris et al. 2000 1. Eu is often thought to 
be synthesized in conditions of very high neutron fluxes (the 
r-process) while Ba is most easily made in conditions of much 
lower neutron fluxes (the s-process). The decline of Ba/Eu may 
show the relative dominance of r-process sites in the early his- 
tory of the Galaxy. 

The first Extremely Metal Poor (EMP, Beers & 
Christlieb I2005i star found to be enriched in n-capture 

elements was the giant HD 115444 with [Fe/H] 3 (Griffin 

et al. 119821 . Westin et al. (2000 1 found marginal overabun- 
dances of Ba and other i-elements relative to iron in this star as 



1 We use the standard notations A/B = N a /Nb and [A/B] = 
log(A' A /A' B ) 1: - \og(N A /N B )o where N x are number densities. 
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compared with the Sun. This still makes the star rich in these 
elements relative to "normal" Population n stars of the same 
metallicity. However, it has a high abundance of europium, 
[Eu/Fe] = 0.85, even compared to the Sun. 

The heavy elements of HD 115444 show good agreement 
in relative abundances with scaled solar r-process abundances, 
as does the EMP giant CS 22892-052 which is even more r- 
element rich (Mc William et al. 1995 Sneden et al. I2003a> . 
This agreement suggests that a similar process is responsible 
for generating the heavy r-elements, both for the most metal- 
poor stars and for the Sun. However, several studies have found 
a deviation from the scaled solar r-process abundance pattern 
for the light r-process elements (e.g. Sneden et al. 2000 2003a 
Cowan et al. 2002 Aoki et al. 2005 1. Also, similar stars like 
CS 31082-001 (Hill et al. IMEl and CS 30306-132 (Honda 
et al. 2004a I seem to have had an "actinide boost", which 
has selectively affected the abundances of the most heavy r- 
elements. These variations in behaviour of the r-elements sug- 
gest that multiple astrophysical sites for the r-process may ex- 
ist. 

The discovery of CS 22892-052 was a result of the compi- 
lation of metal-poor HK survey objects of Beers et al. Q992). 
Similarly, the large programme of Cayrel et al. (2004i at ESO, 
"First Stars", was founded on that survey. In the First Stars 
project, high-resolution studies of the very metal-poor giant 
CS 31082-001 revealed its high Eu abundance. It was possible 
to identify a U n line in this star, and this opened up the possibil- 
ity to use the U/Th ratio as a new chronometer. More recently, 
Honda et al. ( 2004a 2004b I reported on two new Eu-rich stars, 
again taken from the HK survey, namely CS 30306-132 with 
[Eu/Fe] = +0.85 and CS 22183-031 with [Eu/Fe] = +1.2. 
According to the classification suggested by Beers & Christlieb 
(2005 1 the latter of these is an r-II star (having [Eu/Fe] > 
+ 1.0 and [Ba/Eu] < 0.0) while the former is classified as r-I 
(+0.3 < [Eu/Fe] < +1.0 and [Ba/Eu] < 0.0). Another interest- 
ing object which is very overabundant in n-capture elements, 
the carbon-enhanced metal-poor (CEMP) star with s-element 
enhancement CS 31062-050, was recently studied by Johnson 
& Bolte (2004). Interestingly, they concluded that the abun- 
dance profile of the star, in spite of its high Eu abundance, sug- 
gests that the r-process was not responsible in this case. These 
indications will be further discussed in Sect.[6] 

Although the contribution of the s-process to the gas in 
the very early Galaxy may in general be tiny, as one would 
expect due to the small fraction of iron to start the nucle- 
osynthesis from, some stars with [Fe/H] as low as about -3 
have been found to have clear s-element signatures (Johnson 
& Bolte l2002l Simmerer et al. 120041 Sivarani et al. l2004l The 
implications of the presence of these stars are still not fully 
understood. It was predicted by Gallino et al. {T998 1 (see also 
Goriely & Mowlavi 2000 1, that the fraction of the heaviest s- 
elements should be enhanced in metal-poor environments with 
respect to the lighter s-elements, a result of the scarcity of seed 
nuclei with respect to neutrons. This prediction was qualita- 
tively confirmed by Aoki et al. (|2000 1 who made the first dis- 
covery of Pb in a s-element rich very metal-poor star, LP625- 
44. Somewhat more iron-rich "lead stars" were also discov- 
ered by Van Eck et al. (2001) (see also Sivarani et al.|2004). 



We also note that the abundances of the "lighter" n-capture 
elements (Sr, Y, Zr) do not scale very well with those of the 
heavier elements (Ba etc, Z > 56). This may be due to the 
fact that the lighter s-elements (the so-called "weak compo- 
nent", see Prantzos et al. 1990) received important contribu- 
tions from seed nuclei reacting with neutrons created in the 
22 Ne(a',n) 25 Mg reaction in massive AGB stars, while the heavy 
elements were essentially contributed by neutrons from the 
13 C(a,n) 16 reaction in less massive stars. 

Another class of stars was discovered by Barbuy 
et al. ( flW7t and Hill et al. ( ISUUOl Two carbon-rich 
Population n giants, the CEMP stars CS 22948-027 and 
CS 29497-034 from the HK survey, were found not only to be 
rich in the commonly enhanced ^-elements such as Sr, Y, Ba 
and La, but also the r-element Eu was significantly enhanced. 
Preston & Sneden (2001 ) found another star of this class in the 
HK survey (CS 22898-027), Aoki et al. J2UU21 added two more 
(CS 29526-1 10 and CS 31062-012) and Siva rani et al. ( 120041 
another (CS 29497-030) (see also Ivans et al. 12005V The stars 
of this type presently known are listed in Table|8] We shall de- 
note these stars r+s stars, and we postpone a detailed discus- 
sion of the definition of this class of stars to Sect.[6] 

The origin of the abundance peculiarities of the r+s stars is 
not clear, and many scenarios have been proposed (see Sect.|6}. 
However, this is not the only reason for studying them further. 
In fact, the astrophysical sites of the r-process are still unclear. 
The precise sites of the s-process among the most metal-poor 
stars are also still debatable, as well as the fraction of nuclei 
provided by neutrons from the 13 C(a,n) 16 reaction as com- 
pared with the 22 Ne(ff,n) 25 Mg neutron source. One might hope 
that a clarification of the origin of the r+s stars may shed some 
light on the general questions concerning the sites of the r- and 
^-processes. 

The Hamburg/ESO objective-prism survey for bright 
quasars (HES; Wisotzki et al. 2000) has been used for a search 
for metal-poor stars (Christlieb 2003 (. One of these was the 
r+s star HE 2148-1247, found in the Keck pilot programme 
on EMP stars from the HES (Cohen et al. l2"U0ll . In a sys- 
tematic approach to exploring the n-enhanced HES stars, the 
Hamburg/ESO R-process Enhanced Star survey (HERES) has 
been carried out (Christlieb et al. 2004 hereafter Paper i). For 
several hundred candidate stars from HES, "snapshot" spectra 
at moderate resolving power (R ~ 20 000) and low S /N (~ 50) 
were obtained with the main goal to find r-element enriched 
stars, as disclosed by strong Eu lines. These spectra were anal- 
ysed with an automatic procedure based on MARCS model at- 
mospheres and synthetic spectra (Barklem et al. 2005 , Paper n). 
As a result of the HERES survey, a number of new detections 
of n-capture enriched EMP stars were made. In particular, 8 
new r-II stars and 35 r-I stars were found. 

One of the stars found by the HERES survey is 
HE 0338-3945. This object is a star located close to the main- 
sequence turnoff of an old halo population with an overall 

metallicity of [Fe/H] 2.4, and enriched in both Ba and Eu. 

We decided to obtain spectra of higher quality of this star and 
perform a detailed analysis which is presented here. In Sects.|2] 
and[3]the observations and analysis respectively are described. 
In Sect.|4]the abundance results are presented. In Sect.|5]the 
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Table 1. Basic data and parameters for the star HE 0338-3945. 
(') Burstein & Heiles ( f#82t . ( 2 ) Schlegel et al. ( 1998). 



Quantity 


Value 


Uncertainty 


R A (M000) 


0311391^4 9s 


±0.1s 


Decl. (J2000) 


-39°35'44" 


±1" 


V 


15.333 mag 


±0.007 mag 


B-V 


0.420 mag 


±0.016 mag 


V - R 


235 mas 


±0.011 mag 


V - I 


546 map 

V/ ■ J lU 11141.— 


±0.010 mag 


j 


14 403 map 


+0 032 map 

-i-\J .\J J— JUtl^ 


H 


14.1 14 mag 


±0.035 mag 


K 


14.084 mag 


±0.059 mag 


E(B - V)' 


0.000 mag 




E(B - Vf 


0.013 mag 




T c ft 


6160 K 


±100K 


[Fe/H] 


-2.42 dex 


±0.11 dex 


logg 


4.13dex 


±0.33 dex 




1.13km/s 


±0.22 km/s 


v r 


177.9 km/s 


±0.5 km/s 


Class. 


turn-olf star 





abundances are compared with those of other stars. This re- 
veals HE 0338-3945 to be very similar to another r+s star, 
HE 2148-1247 (Cohen et al. l20l)3l . This result led us to sur- 
vey the known r+s stars in the literature, and this survey, along 
with a comparison with HE 0338-3945 is also presented. In 
Sect.|6]the results for HE 0338-3945 are discussed, as is both 
the classification and formation of r+s stars in general. Finally, 
in Sect.0the conclusions are presented. 

2. Observations and data reduction 

The basic data for the star HE 0338-3945 are presented in 
Tabled It was observed in Service Mode at the VLT Unit 
Telescope 2 (Kueyen) equipped with the spectrograph UVES 
at the European Southern Observatory at Paranal, Chile, dur- 
ing the nights of 1 1 December and 23 to 25 December, 2002. 

To cover as large a spectral range as possible, UVES was 
used in dichroic mode with two different settings: Dichroic 1 
with central wavelengths of 3460 A (3055-3850 A) and 5800 A 
(4795-5745 A and 5855-6755 A) in the blue and red arm of the 
spectrograph respectively, and Dichroic 2 with central wave- 
lengths of 4370 A (3760-4945 A) and 8600 A (6720-8425 A 
and 8710-10510A). The exposure times were in total 6 hours 
for the Dichroic 1 setting and 9.5 hours for the Dichroic 2 set- 
ting, consisting of single exposure times of 30 to 75 minutes. 
The projected slit width was set to 1.2", yielding a resolving 
power of R = A/AA equal to 30, 000 - 40, 000. During the ob- 
servations, the CCD binning was set to 1 x 1, but after reduc- 
tion the spectra were rebinned by a factor of two to increase 
the signal-to-noise ratio {S IN) per pixel. Since the pixel scale 
after rebinning is 0.44"/pixel and 0.32"/pixel in the blue and 
red arm, respectively, the spectra are appropriately sampled. 

The optical photometry data (BVRI) for HE 0338-3945 
were obtained on the night of 24 October 2002 by J. Holmberg 
using the Danish 1.5m telescope with DFOSC on La Silla. 
Details of the reduction and analysis of these data are pro- 



vided by Beers et al. (2005, in preparation). These data were 
supplemented by near-IR JHK from the 2MASS Point Source 
Catalog (Cutri et al. l2003l . 

The reduction was made in the standard way by using 
the IDL-package REDUCE, Piskunov & Valenti J2"0"02l There 
were some problems with discontinuities and other reduction 
artifacts, cosmic ray hits, and noise in the spectrum, so the 
usable spectral ranges were finally 3100-5741 A and 5847- 
8487 A, now given in the rest frame of the star. The S /N in 
the spectra was rather poor at the blue end. At 3700 A it had 
risen to about 40, and from 4600 A it was 70 or higher, peaking 
at 6700 A at 130. The S/N was about 74 in mean for the final 
spectra. 

3. Abundance analysis 

3.1. Fundamental parameters and atmospheric model 

In Paper ii, a temperature of 6162 ± 100K was derived from 
photometry when adopting reddening derived from maps of 
Schlegel et al. dl998l >. A subsequent analysis of the snapshot 
spectrum gave logg = 4.09 and [Fe/H] = -2.41. We note that 
if instead the reddening maps of Burstein & Heiles (1982 1 were 
adopted the temperature would be 70 K lower. 

The effective temperature for analysis of the higher quality 
spectrum was redetermined from HB and H6. The continuum 
rectification and analysis, including error estimation, were per- 
formed following Barklem et al. (2002i. These two lines were 
employed as the continuum rectification method was able to be 
successfully applied. The gravity and metallicity from Paper n 
were adopted during the analysis of these Balmer lines. We 
obtained T e ff = 6160 + 140 K, in perfect agreement with the 
result from photometry. Based on these results, we adopted 
an effective temperature of 6160 + 100K for our analysis. 
Attempts were also made to derive the effective temperature by 
comparing the observed colours with calculated colours of the 
MARCS model atmospheres (Edvardsson et al. 2005, in prepa- 
ration). We found values of T e g from B-V, V - R and V - 1 to 
range between 6100 K and 6350 K. 

A small grid of models in log g and [Fe/H] for a fixed T^g 
was calculated with the updated Uppsala model atmosphere 
code MARCS (Gustafsson et al. in preparation). All models 
used scaled solar abundances with the exception of the alpha- 
elements, which were enhanced by 0.4 dex. A microturbulence 
of 1 km/s was adopted for the line opacities. This grid was used 
for an initial analysis of the spectrum, deriving log g and chem- 
ical abundances as described in Sect. 13.31 The parameters from 
Paper ii were used as the initial guess. As C, N and O abun- 
dances are significantly enhanced in this star, the grid was then 
recomputed enhancing these abundances by amounts suggested 
by the initial analysis. The changes in these abundances led 
to a warming of the model atmosphere by about 25 K in the 
typical line forming regions. A final analysis of the spectrum 
was then made as described in Sect. 13.31 giving final model 
parameters of T eB = 6160 ± 100 K, logg = 4.13 ± 0.33 dex, 
and [Fe/H] = -2.42 ±0.11. The final microturbulence £ t was 
1.13 + 0.22 km/s. In Fig.[2we show that adopted model is con- 
sistent with excitation and ionisation equilibrium for Fe and 
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Fig. 1. Plots showing trends of abundances found from individ- 
ual Fe and Ti lines with excitation potential (top), wavelength 
(middle) and equivalent width (bottom). Circles and crosses in- 
dicate lines from the neutral and singly ionised species respec- 
tively, and the dashed lines indicate the determined abundances 
from the simultaneous best fit to all lines of the element in ques- 
tion. 



Ti, and that there are no trends in Fe and Ti abundances with 
line strength or wavelength. We compared these stellar param- 
eters with a [Fe/H] = -2.5 and 12 Gyr isochrone from Kim 
et al. <2002> . and the result lies near the turnoff but between 
the positions of the main-sequence and the sub-giant stage. We 
note that log g is expected to be underestimated since overioni- 
sation of Fe i will lead to an underestimate of the gravity in our 
LTE analysis (estimating the pressure in the atmosphere to be 
lower than it actually is), and thus this star is most likely closer 
to the main-sequence than the sub-giant branch. 



3.2. Line selection and data 



The lines used in the abundance analysis were selected through 
a thorough search and selection process. All lines in the ob- 
served wavelength regions were extracted from the Vienna 
Atomic Line Data Base (VALD, Kupka et al. U999> . Predicted 
equivalent widths were then calculated based on a MARCS 
model atmosphere with the stellar parameters and abundances 
of Paper n. For n-capture element abundances not derived in 
Paper n we used solar values appropriately scaled to have 
similar enhancements to other n-capture elements. Predicted 
strengths of C2, CH and CN features, in the wavelength ranges 
were also calculated. For elements with many available lines 
of reasonable strength, lines with equivalent widths predicted 
to be smaller than < 0.1 mA were rejected, and the remaining 
lines which were predicted not to be overly blended were then 
scrutinised in the UVES spectra, checking that the line is ob- 
served and is indeed unblended. For elements of interest with 
very few lines, (e.g. Ag, Tb) all possible lines were scrutinised. 
The lines from the VALD search were complemented with lines 
of Paper n, and Sneden et al. ( 119961 l27)03a> . The preliminary 
line list was then applied to an initial analysis of the spectrum 
(see Sect. 3.3), and following visual inspection of the results 
each line was either adopted in the final line list or rejected. 
Lines were rejected if the predicted line from the derived abun- 
dance did not agree well with the observed line, with respect 
to the majority lines of the element. This generally indicates 
a blend or erroneous log gf. A large number of lines were re- 
jected in this way. Most blends were found to be due to CH and 
Nd. 



The final line list comprises 621 lines for measuring abun- 
dances and 650 blending lines. Blending lines are modelled 
in the calculated synthetic spectrum using abundances derived 
from other lines, but are not themselves used in determining the 
abundances. Atomic data has been compiled for the selected 
lines. Where several sources of atomic data were available, the 
data were generally extracted in the following order of prior- 
ity: Paper 11, Sneden et al. {2()03ai, Sneden et al. ( 1996), labo- 
ratory data from VALD, Jonsell et al. ( 2005 1, theoretical data 
from VALD. In a few cases, we were aware of recent labora- 
tory /-values which were preferred, e.g. in the case of Nd. In 
general this reflects our preference for well checked data, and 
for laboratory data over astrophysical or theoretical data. The 
original sources of the data were traced where possible, and 
are presented with the final line list in Table 2, which is only 
available electronically. The wavelength A, excitation potential 
X, logg/-value, collisional broadening due to neutral hydro- 
gen, Stark broadening, and radiation damping constant, y ra d, 
are tabulated in Table 2, along with the spectral windows used 
in the abundance analysis (see the description of the method in 
Paper 11). When applicable and available, data for the isotopic 
splitting and hyperfine structure (hfs) are given. Blending lines 
considered are also included in the list. In the Appendix we 
briefly comment on line selection issues and data sources for 
each element. 
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3.3. Analysis method 

The spectrum has been analysed using the automated spec- 
trum analysis code based on SME (Valenti & Piskunov 1996 1 
which is described in detail in Paper n. Here we summarise 
the most relevant aspects of the method. The synthetic spec- 
trum was calculated assuming LTE, a ID plane-parallel geom- 
etry for the atmosphere, and Doppler broadening was modelled 
through the classical microturbulence and macroturbulence pa- 
rameters. Model parameters were optimised to minimise the^f 2 
statistic comparing the synthetic and observed spectra. The ef- 
fective temperature was derived independently as described in 
Sect. l3.U and held fixed throughout. The remaining atmospheric 
parameters (log g, [Fe/H], £ t , v macro ) were first determined from 
an analysis of Fe and Ti lines. These parameters were then 
adopted, and abundances for all elements determined from ap- 
propriate lines considering also blends in the region. As a large 
number of lines regarded as blends were included in this analy- 
sis, a second iteration was performed to ensure that these lines 
were modelled with appropriate abundances. In some instances 
where lines were undetectable, we calculated upper limits to 
the abundances. We computed 3cr upper limits using the same 
code by finding the abundance necessary to produce a line with 
an equivalent width three times the lcr measurement error in 
the equivalent width due to noise that a weak line at the rele- 
vant location in the spectrum would have. 

One significant upgrade of the code was performed, namely 
the partition functions from the MOOG code (Sneden 1973) 
have been adopted, which have been significantly updated and 
corrected in 2002 (see http://verdi.as.utexas.edu/moog.htmli. 
Some significant corrections to the previously used polynomial 
fits to the partition functions of Irwin ( 1981 1 have been found, 
e.g. Tbn (Lawler et al. l2001bi . 

Error estimates in all quantities are computed as detailed 
in Paper n. The propagation of errors from stellar parameters, 
oscillator strengths, observational error and continuum place- 
ment are considered in detail. In Paper n we also estimated 
contributions to the error from modelling uncertainties, such 
as the assumptions of LTE and ID plane-parallel geometry, as 
these may contribute to star-to-star scatter in a sample of stars. 
However, in this case we are examining a single star and this 
error has not been included, i.e. <x e (model) = 0. 

As in Paper n, we make a distinction between absolute and 
relative error estimates. The absolute error, <x abs , estimates the 
total uncertainty in a given quantity. The relative error, cr rel , 
is an error where mainly the uncertainty in the line /-values 
(°"iogg/) nas been neglected. Thus, it is appropriate for star-to- 
star comparison where the same lines and /-values are used, 
as in our homogeneous analysis with HE 2148-1247 (see 
Sect.EHl. 

Since the procedure globally fits the spectral lines of a given 
element, an average uncertainty for each element <r\ ogg f is re- 
quired when computing <x abs . As a number of elements studied 
here were not included in Paper n, in Table[3]we provide an es- 
timated average uncertainty for each element with reference to 
the original literature, which are adopted for o-\ og gf. In the cases 
of Ho and Hf, we could find no error estimate for the oscillator 
strengths and adopted 0.2 dex as an estimate. 
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Table 3. Assigned average values of CHogg/ for each element. 



Element 


""logs/ 


Element 




Element 


""log,?/ 


C 


0.10 


Co 


0.10 


Eu 


0.03 


N 


0.10 


Ni 


0.03 


Gd 


0.05 


O 


0.02 


Cu 


0.04 


Tb 


0.04 


Na 


0.02 


Sr 


0.10 


Dy 


0.04 


Mg 


0.07 


Y 


0.03 


Ho 


0.20 


Al 


0.11 


Zr 


0.03 


El 


0.12 


Ca 


0.11 


Ag 


0.04 


Tm 


0.02 


Sc 


0.04 


Ba 


0.03 


Yb 


0.05 


Ti 


0.05 


La 


0.03 


Lu 


0.04 


V 


0.05 


Ce 


0.10 


Hf 


0.20 


Cr 


0.05 


Pr 


0.04 


Pb 


0.05 


Ma 


0.06 


Nd 


0.03 


Th 


0.02 


Fe 


0.03 


Sm 


0.05 


U 


0.05 



4. Results 

4.1. Abundance results 

We have derived abundances for 33 elements and upper lim- 
its for an additional 6 elements. The results of the abundance 
analysis are given in Table0] and the results for the n-capture 
elements are displayed in Fig. [2] In the table we present the 
derived ID LTE abundances together with the estimated rel- 
ative and absolute errors (see Sect. l3.3> . In calculating the er- 
ror estimates, the sensitivity of each abundance to a change in 
each stellar parameter is calculated, and the results are shown 
in Table|5] 

As has been described above, this analysis is based on the 
assumptions of ID plane-parallel geometry and LTE. Such a 
ID LTE analysis is still the most usual way to estimate the 
chemical composition of stars. However, in recent decades, 
more and more studies of the effects of deviations from these 
traditional assumptions have been carried out, and these stud- 
ies have shown that ID LTE analysis may yield abundances far 
from the "true" values. This is particularly true for metal-poor 
stars, which according to 3D calculations are on average much 
cooler in the line forming layers than predicted by ID model 
atmospheres (Asplund 2004 1. The abundance effects may be as 
large as -0.5 dex for minority species and low excitation lines, 
and even more for molecular lines. Also severe NLTE effects 
may arise due to the steeper temperature structures and temper- 
ature inhomogeneities (Asplund 2005 1. 

While this analysis is a ID LTE analysis, and the ID LTE 
abundances will be discussed in the paper, it is of interest to ex- 
amine the estimated effects of 3D and NLTE, and these should 
be born in mind in interpretation of the abundances. For the 
purposes of our following discussion, we define the estimated 
effects: 

A = log 6nlte,3d — log €lte,id> (!) 

AnLTE = log fNLTE.lD - log eLTE.lD, (2) 

A 3D = loge LT E,3D-logeLTE,iD- (3) 

In the absence of complete line formation calculations in 3D- 
NLTE, of which there are presently very few, the effects for 
the 1D-NLTE and 3D-LTE cases are discussed. Note that A is 
generally different from the sum of Anlte and A3D. In Table[6] 
we present estimates of the 3D and NLTE effects for this star 
and the spectral features employed based on results from the 
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Fig. 2. The abundance pattern of HE 0338-3945 (full circles and downward arrows) compared with the solar r-process pattern 
(solid line) scaled to the Eu abundance of HE 0338-3945, and the solar s-process pattern (dotted line) scaled to the Ba abundance 
of the star. The estimated absolute error bars are shown. The r- and ^-process fractions are from Arlandini et al. ( 1999 1, except 
for Th and U which are from Burris et al. (2000 1. Note, in the case of upper limits (downward arrows) the plotted horizontal bar 
marks the upper limit plus the lcr absolute error estimate in the upper limit. 



literature. The values are in a number of cases uncertain. In 
particular the estimates of NLTE effects are often plagued by 
uncertainties in collision rates. Note that a range of 3D effects 
are often quoted as the effects on individual lines may vary. 

Below, we discuss the results and relevant issues in the 
analysis for each element in turn. Elements not discussed in- 
dividually have a significant number of unblended lines with 
adequate S /N for a clear detection were available, the analysed 
lines were all well fit by a single abundance using our adopted 
model, no NLTE or 3D corrections were available, and the 
abundance was unremarkable with respect to "normal" metal- 
poor stars (see Sect. 15. 2> . Note that line selection and atomic 
data for each element are discussed in the Appendix. 

Carbon: HE 0338-3945 is carbon-rich and the two strong 
bands and two lines of CH used in the abundance analysis were 
well fitted. The abundance is, however, uncertain, due to the 
strong temperature sensitivity of the CH molecule formation 
(see Table[5}- In addition, preliminary calculations suggest that 



A3D may amount to approximately -0.5 to -0.3 dex (Asplund 
& Garcia Perez l2U0Tl Asplund 20©. 

The isotopic ratio 12 C/ 13 C was determined from reasonably 
well isolated 13 CH features between 4210 and 4225 A. The 
comparisons of observed and synthetic spectra are shown in 
Fig. [5] where the total derived C abundance derived is adopted. 
We find 12 C/ 13 C~ 10. 

Nitrogen: The N abundance was derived based on bands of 
the CN molecule. The abundance is uncertain, first because it 
is both highly temperature sensitive (see Table [5J and depen- 
dent on the C abundance, and secondly because it is probably 
subject to 3D effects. The value of A3D for CN is expected to be 
lower than for NH, which may be as much as -0.6 to -0.9 dex 
for turn-off stars (Asplund 2005 1. A reasonable estimate of the 
effect may amount to -0.5 dex. Note, the NH An-XZ band at 
3360 A gave an abundance consistent with the CN result, but 
was disregarded in the final analysis due to poor S /N and a 
large degree of blending. 
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Table 4. Derived elemental abundances for HE 0338-3945. For each element X the LTE abundance is presented as log ex and 
[X/Fe], together with their relative, cr rel , and absolute, cr abs , r.m.s. error estimates (see text). N, , gives the number of lines of 
the element used for abundance determination for each element, and gives the number of those features classified as 3cr 
detections. 



Element 




log 6x 


erf 


log fx 


[X/Fe] 


o- rel 

[X/Fe] 


abs 

[X/Fe] 


Nat 


N 3(r 








[dex] 


[dex] 


[dex] 


[dex] 


[dex] 


[dex] 






6 


C 




8.08 


0.09 


0.15 


2.13 


0.06 


0.15 




_l 


7 


N 


(TN> 


6.92 


0.16 


0.20 


1.55 


0.09 


0.17 




_i 


8 


O 


(Ci T^l 


7.66 


0.04 


0.09 


1.40 


0.09 


0.1 1 


2 


1 


11 


Na 


(Na TI 


4.26 


0.03 


0.06 


0.36 


0.07 


0.10 


3 


2 


12 


Mg 


Vivig i; 


5.45 


0.02 


0.09 


0.30 


0.09 


0.14 


7 


5 


13 


Al 


(Al TI 


3.16 


0.09 


0.15 


-0.88 


0.03 


0.13 


1 


1 


20 


Ca 


(Ca Fl 


4.31 


0.05 


0.12 


0.38 


0.07 


0.14 


10 


7 


21 


Sc 


(Sc II) 


1.26 


0.13 


0.16 


0.53 


0.05 


0.09 


5 


2 


22 


Ti 


(Ti I+II) 


2.96 


0.12 


0.14 


0.37 


0.04 


0.08 


24 


1 1 


23 


V 


(V I+II) 


1.76 


0.13 


0.16 


0.19 


0.05 


0.09 


5 


o 


24 


Cr 


(Cr I+II) 


3.12 


0.09 


0.1 1 


-0.12 


0.03 


0.08 




2 


25 


Mn 


(Mn T+TH 


2.47 


0.11 


0.12 


-0.49 


0.04 


0.08 


5 


4 


26 


Fe 


(Fe I+II) 


5.07 


0.10 


0.11 


0.00 


0.03 


0.05 


69 


30 


27 


Co 


(Cn Fl 


2.72 


0.09 


0.14 


0.23 


0.03 


0.11 


4 


3 


28 


Ni 


(N'\ i) 


3.83 


0.10 


0.1 1 


0.01 


0.03 


0.06 




5 


29 


Cu 


(Cu I) 


1.02 


0.12 


0.12 


-0.75 


0.06 


0.08 


2 


o 


38 


Sr 


(Sr I+II) 


1.28 


0.08 


0.16 


0.74 


0.03 


0.14 


5 


2 


39 


Y 


(Y TT1 


0.63 


0.13 


0.16 


0.83 


0.05 


0.08 


1 1 


6 


40 


Zr 


(7x rn 


1.37 


0.13 


0.15 


1.20 


0.04 


0.08 


16 


7 


47 


Ag 


(Ac TI 


< 0.45 


0.15 


0.16 


< 1.94 


0.11 


0.12 


2 


o 


56 


Ba 


(Ba II) 


2.10 


0.06 


0.07 


2.41 


0.05 


0.08 


3 


3 


57 


La 


(La II) 


1.01 


0.13 


0.16 


2.28 


0.04 


0.08 




6 


58 


Cc 


(Ce II) 


1.30 


0.14 


0.18 


2.16 


0.05 


0.12 


15 


g 


59 


Pr 


(Pr TTl 


0.44 


0.16 


0.19 


2.16 


0.07 


0.11 


4 


o 


60 


Nd 


(Nd II) 


1.37 


0.14 


0.16 


2.30 


0.05 


0.08 


34 


15 


62 


Sm 


(Sm II) 


0.82 


0.15 


0.18 


2.25 


0.06 


0.10 


4 


1 


63 


Eu 


(Eu II) 


0.02 


0.14 


0.17 


1.94 


0.05 


0.09 


3 


3 


64 


Gd 


(Gd II) 


0.69 


0.13 


0.16 


2.00 


0.05 


0.09 




o 


65 


Tb 


(Th m 


< 0.13 


0.16 


0.18 


< 2.66 


0.09 


0.11 


1 


o 


66 


Dy 


(Dy II) 


0.55 


0.14 


0.17 


1.84 


0.05 


0.09 


16 


3 


67 


Ho 


(Ho II) 


< 1.16 


0.15 


0.27 


< 3.33 


0.08 


0.23 


1 





68 


Er 


(Er II) 


0.74 


0.14 


0.20 


2.24 


0.05 


0.15 


3 


1 


69 


Tm 


(Tm II) 


0.06 


0.13 


0.15 


2.49 


0.05 


0.08 


4 





70 


Yb 


(Yb II) 


0.77 


0.13 


0.15 


2.12 


0.07 


0.09 


2 


2 


71 


Lu 


(Lu II) 


<0.72 


0.15 


0.17 


<3.09 


0.09 


0.12 


1 





72 


Hf 


(Hf II) 


0.94 


0.17 


0.28 


2.49 


0.10 


0.23 


1 


1 


82 


Pb 


(Pbl) 


2.62 


0.10 


0.11 


3.10 


0.04 


0.08 


2 


1 


90 


Th 


(Th II) 


<0.23 


0.17 


0.19 


< 2.57 


0.08 


0.11 


1 





92 


U 


(U II) 


< -0.11 


0.38 


0.74 


< 2.82 


0.46 


0.80 


1 






(') The detection exceeds 3<r for the molecular bands seen as single features. 



Oxygen: Only two of the three lines of the infrared O i triplet 
at 7773 A were used in the analysis, as the 7775 A line was 
affected by a reduction artefact. The [Oi] line at 6300 A 
could not be detected. A value of A NLTE ~ -0. 1 dex has 
been adopted for stars with similar fundamental parameters 
by Nissen et al. ( I2002> . A3 D is estimated to be approximately 
-0. 1 dex by Asplund ( 120051 . 



Sodium: The derived Na abundance was dominated by the 
Nai D lines at 5890 A, and is determined quite precisely as 
seen from the small estimated errors in Table|4] However, the 



true Na abundance may be lower due to NLTE effects. From 
interpolation in Table 2 of Baumiiller et al. ( 1998) we estimate 
Anlte ~ -0.4 dex. The 3D effect seems not to be severe for 
sodium (Asplund 20041. 



Magnesium: The lines of Mgi were well fit, and the er- 
ror in the abundance estimated to be quite small. However, 
it is affected by both NLTE and 3D effects. According to 
Asplund ( 2005 1 the NLTE effect on the abundance from the 
Mgi line at 5 172 A is around -0.2 dex, while for other Mgi 
lines the effect is approximately +0.2 dex. We have used a total 
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Table 5. The change in the logarithmic abundances due to a 
change of the fundamental parameters by amounts correspond- 
ing approximately to their uncertainties. 



Element 


r cl r + 100 K 


logg + 0.3 dex 


f i +0.2 km/s 




[dex] 


[dex] 


[dex] 


C 


0.15 


-0.09 


-0.01 


N 


0.23 


-0.09 


0.00 





-0.06 


0.10 


0.00 


Na 


0.06 


-0.04 


-0.03 


Mg 


0.06 


-0.06 


-0.01 


Al 


0.1 1 


-0.03 


-0.03 


Ca 


0.06 


-0.02 


0.00 


Sc 


0.05 


0.10 


-0.03 


Ti 


0.05 


0.08 


-0.03 


V 


0.05 


0.09 


0.00 


Cr 


0.09 


-0.01 


-0.04 


Mn 


0.10 


0.01 


-0.01 


Fe 


0.07 


0.03 


-0.02 


Co 


0.09 


0.01 


-0.01 


Ni 


0.09 


0.01 


-0.04 


Cu 


0.10 


0.00 


-0.01 


Sr 


0.08 


-0.02 


-0.09 


Y 


0.05 


0.09 


—0.04 


Zr 


0.06 


0.09 


—0.02 


A 1 

Ag 


0. 1 1 


0.01 


-0.01 


Ba 


0.08 


-0.03 


-0.02 


La 


0.05 


0. 10 


—0.01 


Ce 


0.06 


0.10 


-0.01 


Pr 


0.06 


0.12 


0.00 


Nd 


0.06 


0. 10 


-0.02 


Sni 


U.U / 


n i i 
0. 1 1 


0.00 


En 
till 


n oa 


n i i 

U. 1 1 


o on 
u.uu 


Gd 


0.05 


0.10 


0.00 


Tb 1 


0.07 


0.09 


-0.01 


Dy 


0.07 


0.09 


-0.02 


Ho 1 


0.05 


0.11 


0.00 


Er 


0.06 


0.10 


-0.01 


Tm 


0.05 


0.09 


-0.01 


Yb 


0.07 


0.06 


-0.03 


Lu 1 


0.05 


0.09 


0.00 


Hf 


0.08 


0.09 


-0.06 


Pb 


0.11 


-0.02 


-0.01 


Th 1 


0.09 


0.10 


-0.01 


U 1 


0.35 


-0.80 


-0.01 


(' ) Effects on the upper limits are shown. 



of 7 lines, and so Anlte should be positive. This is supported by 
Gehren et al. (2004i who show results for 7 lines, although not 
exactly the same set of lines as used here. The resulting value 
of Anlte for stars with parameters close to HE 0338-3945 is 
~ +0.1 dex. Asplund (2004i estimate that A 3D may be -0.5 to 
-0.3 dex for the low excitation lines. 

Aluminium: The abundance of Al was derived using only the 
strong resonance line of Al i at 3961 A, and gives a clearly sub- 
solar value relative to Fe as is the norm for metal-poor stars (see 
Sect. 15.21 . Anlte could, according to an interpolation made in 
Table 1 of Baumuller & Gehren (119971 . be as high as +0.5 dex. 
According to Asplund ( 20041 120051 A 3 d may amount to a sim- 
ilar value but of opposite sign for low excitation lines. 

Calcium: The calcium abundance was determined from lines 
of Ca i. A value of Anlte ~ +0.2 dex was adopted from re- 
cent work of Korn & Mashonkina (2006, in preparation), and 
according to Asplund (2004i the A 3 d may amount to -0.3 to 
-0.1 dex. 



Table 6. Estimated effects of NLTE and 3D on abun- 
dances from literature and adapted to the parameters of 
HE 0338-3945. The numbers are in many cases highly 
uncertain; see text for further discussion. References: (a) 
Asplund (120041 (b) Asplund (120051 . (c) Asplund & Garcia 
Perez d2"U0"Tl (d) Baumuller et al. ©9Si», (e) Baumuller & 
Gehren ( fl997l . (f) Collet et al. (120051 . (g) Gehren et al. (120041 
(h) Korn & Mashonkina (2006, in preparation) (i) Nissen 
et al. (120021 . 



Element 


Anlte 
[dex] 


Ref-NLTE 


Aid 
[dex] 


Ref,, D 


6 C (CH) 






-0.5 to -0.3 


a.c 


7 N (CN) 






-0.5 


see text 


8 O 


-0.1 


i 


-0.1 


b 


11 Na 


-0.4 


d 


-0.1 to +0.1 


a 


12 Mg 


+0.1 


g 


-0.5 to -0.3 


a 


13 Al 


+0.5 


e 


-0.5 to -0.3 


a 


20 Ca 


+0.2 


h 


-0.3 to -0.1 


a 


26 Fe (Fe i) 


+0.6 


f 


-0.5 to -0.3 


a 


26 Fe (Fe n) 


+0.2 


b 


-0.1 to +0.1 


a 


38 Sr 


+0.3 


b 


-0.5 to -0.3 


a 


56 Ba 


+0.2 


b 


-0.5 to -0. 1 


a 


63 Eu 


+0.1 


b 








4218.0 4218.5 4219.0 4219.5 

A [A] 

Fig. 3. Comparison of observed and synthetic spectra for 13 CH 
features. The double line shows the observed spectrum. The 
thin lines show synthetic spectra for isotope ratios I2 C/ I3 C=3, 
10, and 30. 
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Scandium: The abundance of scandium was determined using 
5 Sc ii lines. The abundance seems to be high in comparison to 
normal metal-poor stars (see Sect. l5.2l i. The two bluest lines are 
strongest and dominate the fitting, and are thus well fit. The fit 
does, however, seem to overestimate the strength of the three 
weak lines redward of 5000 A. This may suggest that the blue 
lines are affected by unknown or poorly modelled blends. We 
compared the [Sc/Fe] and [C/Fe] abundances from Paper n, and 
found that while the results are uncorrelated for [C/Fe] i= 1.5, 
for [C/Fe] <: 1.5 there is a trend to high [Sc/Fe] with high 
[C/Fe]. Thus, our Sc abundance may be overestimated. 

Vanadium: Five weak lines in the blue were used to derive 
the abundance, the best being only detected at the 1.8<x level. 
However, all lines of both species, V i and V n, are well mod- 
elled and we consider this a reliable detection. 

Manganese: Three lines of Mnn and two of Mm were anal- 
ysed. The S IN was significantly better at the Mm lines than 
at the Mnn lines, and thus the Mm lines dominated the fit. 
However, the derived abundance did not reproduce the Mnn 
lines well, the model prediction being significantly weaker than 
the observed by a consistent amount for all three lines. The dis- 
crepancy may be due to overionisation of Mn i. Such a discrep- 
ancy has also been noticed by Johnson (2002 1. 

Iron: A total of 61 Fei and 8 Fen lines were analysed. In 
Table[6]we give estimates for A^d and Anlte for each species 
separately. As our Fe abundance is dominated by the numerous 
Fe i lines the effects on Fe i should dominate. Note, however, 
that the effects of NLTE are presently highly uncertain due to 
uncertainty in inelastic collision rates due to hydrogen atoms. 
For example, Korn et al. (2003 1 found much smaller NLTE ef- 
fects when astrophysically calibrating these collision rates. As 
noted in Sect. B.ll NLTE effects on Fe i are expected to lead to 
underestimation of log g in LTE, and this would affect all other 
elements. 

Copper: The Cu abundance was derived from two lines in the 
ultraviolet. Both lines were clearly seen and identified, and well 
modelled with a single abundance. However, the low S /N ratio 
in this spectral region meant that the strongest line was detected 
only at the 2.9<x level, and continuum placement was uncertain. 
Despite this, we regard it as a definite detection. 

Strontium, yttrium, zirconium: Sr, Y andZr, are strongly over- 
abundant compared to normal metal-poor stars (see Sect. l5.2> . 
Interestingly, the LTE abundances of these elements lie be- 
low the solar s-process abundance pattern as normalised to Ba 
(see Fig. 0, and curiously enough, at least Y and Zr fit the re- 
process abundance pattern normalised to Eu. Strontium is es- 
timated to have Anlte ~ +0.3 dex for the lines at 4077, 4161 
and 4215 A in metal-poor turn-off stars (Asplund 2005 1. We 
analysed two more weak lines, one with low S IN. The stronger 
lines at higher S /N have high weight in the abundance calcula- 



tion, and therefore this estimate for Anlte is appropriate. A3D 
is approximately -0.3 to -0.5 (Asplund l2004> . 

Silver: The upper limit to the silver abundance is estimated 
from two Ag 1 lines, both blended and of roughly equal intrinsic 
strength. The upper limit is sufficiently low to be inconsistent 
with the scaled solar r-process abundance pattern (see Fig.[2}- 
This may reflect overionisation of Ag 1 (see Sect.|6j. 

Barium: Ba is strongly overabundant compared to normal 
metal-poor stars (see Sect. 15. 2> . This abundance was derived 
using three Ban lines: two weak and the strong resonance line. 
The resonance line at 4554 A has an estimated value of Anlte 
of at least +0.2 dex (Asplund 2005 1. A3D is of the opposite sign, 
and is estimated at -0.5 to -0.1 dex (Asplund 20041. 

Lanthanum: Lanthanum is also strongly overabundant in this 
star. Eight lines of La 11 of varying strength were employed, 
the majority of which were well fit by a single abundance. 
However, two lines, 4920 and 4921 A were not well fit, the 
observed lines being stronger than the predicted lines for the 
derived abundance. We note that for these lines we had incom- 
plete hfs data (see Appendix), and this could possibly lead to 
the line strengths being underestimated. 

Praseodymium, Neodymium, Samarium: Four weak lines 
were employed to derive the abundance of Pr, the best being 
only detected at the 1 .6<x level. However, all 4 lines were well 
modelled, and thus we consider this a reliable detection. A total 
of 34 mostly weak lines were used to derive the Nd abundance, 
only slightly less than half of the lines being detected at the 3cr 
level. For Sm, four weak lines of Sm 11 were used. 

Europium: The abundance of Eu is enhanced compared to 
typical metal-poor stars (see Sect. l5.2> . The abundance may 
be even higher, as Anlte is probably greater than +0.1 
(Asplund 2005 1. No estimate of A3D is available, but see the 
discussion regarding abundance ratios at the end of this sec- 
tion. 

Gadolinium: Eight weak lines are used, the best being only a 
1.8cr detection. However, all eight lines are consistently mod- 
elled, and we class this as a reliable detection. 

Terbium: No line was detected, and thus an upper limit was 
derived from single line of Tb 11. 

Holmium: No lines are detected, and thus a 3<x upper limit was 
derived. Due to strong blending at the strongest lines of Ho, we 
were forced to use a rather weak line and thus the derived upper 
limit is rather high and a weak constraint. 

Erbium, Thulium: The abundance of Er was based on three 
weak Em lines. For Tm, 4 weak lines were employed, the best 
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being only detected at the 1.6cr level. However, the lines of 
these elements were well modelled, and thus we consider these 
as reliable detections. Both elements have abundances signif- 
icantly above the scaled solar s- and r-process patterns (see 
Fig-E}. 

Lutetium: No lines were detected, and thus a 3<x upper limit 
was derived. The employed line region is somewhat blended, 
but we regard the derived upper limit as reliable. 




10 20 30 40 50 60 
atomic number Z 



Hafnium: Only the resonance line in the ultraviolet was avail- Fig. 4. Difference between the abundances for HE 0338-3945 
able, but a definite detection was made. from this work and those from Paper n. 



Lead: The Pb abundance is significantly above the scaled s- 
process abundance pattern in Fig.|2] Pb has an ionization en- 
ergy of 7.4 eV, and the two lines observed are both Pbi lines. 
Here, overionisation is a possibility, which would lead to an un- 
derestimated Pb abundance and e.g. Pb/Ba ratio, meaning that 
this ratio, already remarkably high in the LTE analysis, may be 
even higher. 

Thorium: The 3<x upper limit was based on a single Thu line 
at 4019.129 A. The region contains a number of blends, but a 
reliable upper limit was possible to derive. 

Uranium: No lines are detected, so a 3<x upper limit has been 
derived. The region is significantly blended, and the reliability 
of the derived upper limit is dependent on the correct modelling 
of these blends. 

Some specific abundance ratios: We shall finally comment 
on the accuracy of some important abundance ratios, to be fur- 
ther discussed below. The similarity of Sr, Y, Zr and Ba in terms 
of ionisation energies and atomic structure seems to suggest 
that the effects of NLTE should be limited on Ba/Sr, Ba/Y and 
Ba/Zr. Also, the most important lines in the analysis have low 
excitation energies - the Ba abundance determination is dom- 
inated by the wings of the 4554 A resonance line, while the 
other two lines with 2.5 and 2.7 eV excitation energies are less 
significant. Thus, the temperature sensitivity of the important 
line strengths are roughly similar, which is reflected in simi- 
lar temperature sensitivities of the abundances in Table|5] This 
should lead to similar sensitivities to 3D effects. Similarly, the 
effects of NLTE and 3D on abundance ratios of these elements 
relative to Eu may be expected to be small. 

4.2. Binarity 

Radial velocities were measured from a series of Fe i lines with 
accurate laboratory wavelengths between 4000 and 5000 A, us- 
ing the snapshot spectrum observed on 15 October 2002, and 
the observations from December 2002 which comprise the high 
quality spectrum analysed in this work. No difference in ra- 
dial velocity shift was observed between the different expo- 
sures taken in December 2002. The barycentric radial velocities 



measured from the snapshot spectrum and the co-added high 
quality spectrum were 177.7 + 0.8km/s and 177.9 + 0.5km/s 
respectively, where the quoted errors are the standard devia- 
tions of the results from different lines; errors due to the wave- 
length calibration uncertainties are not included. We therefore 
presently have no evidence for the binarity of HE 0338-3945. 
A recent study by Lucatello et al. (|2005 1 suggests all CEMP 
stars are multiple systems. 

5. Comparisons 

5.1. Comparison with the results of Paper u 

As has been mentioned, a snapshot spectrum of this star was 
analysed in Paper n. The snapshot spectrum of HE 0338-3945, 
compared to the spectrum employed in this work, is of rela- 
tively low quality with S /N ~ 47 and R ~ 20000 and has sig- 
nificantly less wavelength coverage, A = 3760-4980 A. The 
snapshot spectrum was reduced by the UVES pipeline, while 
the spectrum analysed here was reduced as described in Sect.|2] 
However, the spectra have been analysed in a similar manner, 
and it is of interest to check if the results are consistent. We 
have already seen in Sect.[3] that the results for the stellar pa- 
rameters are in good agreement. Figure compares the abun- 
dances where available from Paper n, and we see that the results 
are in good agreement. 

5.2. Comparison with normal metal-poor stars of 
Paperu 

To place HE 0338-3945 in context, it is of interest to compare 
this star with "normal" Population n stars of similar metallicity. 
We extracted all stars from Paper ii with -2.6 < [Fe/H] < -2.2. 
It was required that the stars should have [Eu/Fe] < 1 and 
[C/Fe] < 1 to be considered as "normal". The mean value of 
[X/Fe] was then computed for the sample where more than five 
measurements were available. 

In Fig. [5] we compare the abundances for HE 0338-3945 
with the described mean abundances. For reference, the de- 
rived mean abundances and standard deviations are provided 
in Table0 In Fig.|5] we clearly see the enhancements of C 
and the «-capture elements relative to the normal stars. For the 
most part the iron-peak element abundances for HE 0338-3945 
appear to be normal; however, Sc appeai-s to be significantly 
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Fig. 5. Difference between abundances for HE 0338-3945 and 
those for normal Population n stars of similar metallicity from 
Paper ii, derived as described in the text. The vertical bars show 
the standard deviation among the normal Population n stars. 

Table 7. Comparison of the mean abundances relative to Fe, 
for normal Population n stars of Paper n in the metallicity range 
-2.6 < [Fe/H] < -2.2, with those for HE 0338-3945. For each 
element we report the number of measurements N available 
from Paper n, the mean value of [X/Fe], and the standard devia- 
tion of [X/Fe] for the sample, the abundance for HE 0338-3945 
[X/Fe]*, and the difference [X/Fe]* - ([X/Fe]). 



Element 


N 


<[X/Fe]> 


o-tfX/Fe]) 


[X/Fe], 


A[X/Fe] 


C 


58 


0.18 


0.30 


2.13 


1.95 


Mg 


58 


0.26 


0.08 


0.30 


0.04 


Al 


55 


-0.86 


0.23 


-0.88 


-0.02 


Ca 


58 


0.27 


0.10 


0.38 


0.11 


Sc 


58 


0.08 


0.13 


0.53 


0.45 


Ti 


58 


0.24 


0.12 


0.37 


0.13 


V 


18 


0.00 


0.18 


0.19 


0.19 


Cr 


58 


-0.19 


0.14 


-0.12 


0.07 


Mn 


58 


-0.47 


0.15 


-0.49 


-0.02 


Co 


53 


0.13 


0.10 


0.23 


0.10 


Ni 


58 


-0.06 


0.14 


0.01 


0.07 


Sr 


58 


-0.10 


0.29 


0.74 


0.84 


Y 


35 


-0.20 


0.16 


0.83 


1.03 


Zr 


16 


0.27 


0.18 


1.20 


0.93 


Ba 


55 


-0.30 


0.42 


2.41 


2.71 


La 


6 


-0.06 


0.24 


2.28 


2.34 


Nd 


10 


0.24 


0.52 


2.30 


2.06 


Eu 


21 


0.39 


0.29 


1.94 


1.55 



enhanced compared to the normal stars. The good agreement 
found in Sect. 15. II with results for this star from Paper n means 
that this cannot be due to a difference in the analysis, e.g. 
use of different lines or atomic data. However, as discussed in 
Sect. 14. li the Sc abundance may be spurious. 

5.3. Comparisons with the star HE 2148-1247 

Soon after examining our results for HE 0338-3945, it 
became apparent that this star was very similar to the 
star HE 2148-1247 studied by Cohen et al. J20Q31 . The 
Keck/HIRES spectrum of this object used by Cohen et al. 
was kindly provided to us by Judith Cohen. Order merg- 
ing and continuum placement was performed by us; all other 
reduction steps had been performed. Thus, using our auto- 
mated spectrum analysis technique we were easily able to 
perform an analysis of this spectrum homogeneous with our 
analysis of HE 0338-3945. Note, as the spectral coverage 



Fig. 6. Difference between abundances for HE 0338-3945 (cut 
spectrum) and HE 2148-1247. The dashed line shows the 
mean difference. 

of the HE 2148-1247 spectrum was not as extensive as for 
our HE 0338-3945 spectrum, we produced a spectrum of 
HE 0338-3945 with identical wavelength coverage to the 
HE 2148-1247 spectrum, which we refer to as the "cut" spec- 
trum of HE 0338-3945. Thus, an analysis of the cut spectrum 
permits a more strictly homogeneous analysis as exactly the 
same spectral lines are employed. 

For our analysis of HE 2148-1247, we adopted the effec- 
tive temperature from Cohen et al., T e g = 6380 K. The abun- 
dance results from our analysis of HE 2148-1247 and the cut 
spectrum of HE 0338-3945 are compared in Fig.[6] The abun- 
dance patterns of the two stars are seen to be very similar. An 
offset of -0.2 dex between the abundances is seen, but this 
could easily be due to errors in T e s which were not derived 
consistently. The scatter is consistent with the uncertainties. 
From this we conclude that the stars have identical abundance 
patterns for the elements considered, within the errors in the 
analysis. 

5.4. Comparisons with other similar stars 

In addition to HE 0338-3945 and HE 2148-1247 there are 
several other stars that have recently been found to show both 
enhanced r- and s-elements. In total 17 such stars are known to 
us, and these are listed in Table|8] 

We have investigated the similarities within the r- and s- 
element enhanced group, and the differences relative to other 
metal-poor stars. The r- and s-enhanced stars were compared 
with other stars of [Fe/H]< -1.0 investigated in Paper n, by 
Burris et al. (2000 1, and Fulbright (2000 1, and some other 
groups listed in Table|S] When several sources were available, 
data were adopted from the source considered most reliable or 
the mean of recent data was taken. 

The r- and s-enhanced stars HE 1046-1352 and 
HE 1405-0822 were analysed in the work of Paper n, but un- 
published due to the significant pollution of the spectrum by 
molecular carbon features making the analysis method uncer- 
tain. The presented data for these stars must therefore be re- 
garded as uncertain and preliminary. However, higher quality 
data of the star HE 1405-0822 is considered in the HERES pa- 
per Sivarani et al. (in preparation). It must also be observed that 
the error bars were quite large for some stars in several investi- 
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Table 8. Data for stars of different classes. For stars r-I stars and stars not n-enhanced, see Paper n, Bums et al. (2000), and 
Fulbright (2000), although data are presented for some outliers of interest from Burris et al. When several investigations exist, 
either the most reliable source or the mean of recent data has been taken. The mean of abundance ratios [X/Y], Y^Fe, was 
calculated from the mean of the abundances relative to Fe, and not from individual measurements of [X/Y]. (:) indicates an 
insecure measurement, and such cases are not considered when calculating the means. References: (1) This work, (2) Aoki 
et al. ( 120011 . (3) Aoki et al. (ED02), (4) Barbuy et al. (2005 1, (5) Paper n (6) Unpublished preliminary result obtained in Paper n, 
(7) B urris et a l. (120001 . (8) Coh en et a l. (120031 . (9) Hill et al. d20U0l , (10) Hill et al. |20Q1 a nd Ple z et al. ( 120041 . (11) Hond a 
et al. d2004a> . (12) I vans e t al. (120051 . (13) Johnso n & Bo lte (120021 . (14) Joh nson & Bolte d2004l . (15) Jonsell et al. J20051 . 
(16) Lucatello et al. (120031 . (17) Preston & Sneden HM\\ . (18) Sivarani et al. d2004l . (19) Sneden et al. ( I2003a> (20) Van Eck 
et al. ( 120031 . (21) Zacs et al. JW981 . (22) mean of presented data, (23) mean of (4) and (17). 
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0.26 


-0.13 






(7) 


BD+11°2998 


5425 


2.30 


-1.17 




-0.20 


-0.23 


0.17 


0.14 


-0.17 




0.06 




0.08 


-0.23 






(71 


BD+17°3248 


5250 


2.30 


-2.02 




0.55 


0.06 


0.50 


0.97 


0.60 




0.96 




0.01 


-0.36 






(7) 
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gations, and that errors in stellar parameters may be the reason 
for stars with seemingly peculiar abundance patterns. 

The plot in Fig.[7]shows a clear separation of different kinds 
of stars. The criteria of the classes in Table|5]in Sect.[6] were 
defined to achieve this. The stars were classified as r-I, r-II, s- 
element enhanced, r+s, i.e. both r- and ^-element enhanced, 
and/or Pb enhanced. Stars with no significant enhancement, 
and thus not falling under any of these criteria, are in this paper 
called "normal" metal-poor stars. 

The r+s stars are rich in both r- and s-elements, resid- 
ing in the upper right corner of Fig.0 All of these stars have 
[Ba/Eu] > 0.0. This is also true for the Pb stars with mea- 
surements of both Ba and Eu found in the literature. These Pb 
stars are also then classified as r+s stars. The ^-enhanced stars 
also have [Ba/Eu] > 0.0, and may be located along the same 
sequence in the figure as the r+s stars. The r-II and r-I stellar 
classes both have [Ba/Eu] < 0.0. The r-II stars have in gen- 
eral a typically lower content of [Eu/Fe] than r+s stars, and in 
Fig-d seem to separate clearly from the r+s stars. The gap in 
between the groups might be due to poor statistics; this sep- 
aration was also noted by Cohen et al. ( 2003 ) (however, with 
some overlap). The normal stars, i.e. stars with no n-capture 
element enhancement, have predominantly [Ba/Eu] < 0.0, al- 
though positive ratios may also be found. The same pattern was 
also found by Cohen et al. ( 2003 Fig. 13) and Johnson & Bolte 
(2004 Fig. 5), although the nomenclature is different in these 
studies (see Sect.Q. 

The plots in Fig.|S] show the ratio of the different heavy 
n-capture s-elements (hs = Ba, La, Ce) compared to the r- 
element Eu for stars of different classes. As is also seen in Fig.0 
there is a clear separation between on one hand r+s and s stars, 
and on the other hand r-II, r-I and normal stars. The different 
classes of stars seem to follow a single well-defined sequence in 
the [La/Eu] vs. [Ce/Eu] plot, but less so in the other abundance 
diagrams. [Ba/Eu] shows the greatest variation amplitude in 
the diagrams, while La and Ce seem vary less and follow each 
other. According to Paper n this may well be an observational 
effect; as seen in Table[5]the abundances of Ba and of the other 
elements react differently to changes in effective temperature 
and gravity which is essentially a consequence of the signifi- 
cance of the strong damping wings of the dominating Ba n line 
/14554. Therefore, stellar-parameter errors cancel for the ratios 
[La/Eu] and [Ce/Eu], but not for [Ba/Eu]. The Pb enhanced 
r+s star CS 29526-1 10 has a rather low abundance of La (Aoki 
et al. 12002b . which makes the ratio [La/Eu] negative. This is 
astonishing since the other ns-elements (Ba, Ce) are consider- 
ably enhanced. The authors note that the effective temperature 
and gravity are rather uncertain for this star. 

There are some outliers among the normal stars in the upper 
panel of Fig.[8] The parameters of these stars may be viewed 
in Table[8] These stars reside in or quite near the same region 
as the r+s stars, but are neither r- nor s-enhanced. The star 
BD+17°3248 is rather close to being regarded as an r+s star, 
as it has abundances [Ba/Fe] = 0.97, [Eu/Fe] = 0.96, and 
[Ba/Eu] = 0.01. However, in Fig.0it lies closer to the nor- 
mal stars. We note that Cowan et al. 2002 find a much lower 
Ba abundance for this star, which would place it firmly as a r-I 
star. 
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Fig. 7. A plot comparing the abundance of the s-element Ba to 
the r- element Eu for stars of different classes (see definitions 
in Table|5}. The r+s stars, enhanced in both r- and s-elements, 
are shown as open circles, and the r+s stars also known to be 
enhanced in Pb are shown as circles with asterisks inside. The 
stars enriched only in r-elements are marked as small (r-I) and 
large (r-II) open squares, and stars with only ^-enhancement 
are shown as diamonds. Stars not enhanced in either r- or s- 
elements, in this article called normal stars, are marked by 
small filled circles. The dotted lines mark one-to-one relations. 
Data are taken from Tableland for normal and r-I stars from 
Paper n, Burris et al. (2000), and Fulbright d^UoDV 

The histograms in Fig.|9]display the abundance of the heavy 
s-element Ba compared to the light ^-element (Is) Sr for the 
different stellar classes. In spite of the small number statistics, 
there are clear separations in mean ratios for different classes. 
Normal stars have a slightly negative mean [Ba/Sr], although 
the spread is wide. The mean r-I star has slightly more Ba than 
Sr, and for r-II stars Ba dominates. The s stars in the diagram 
are very few, but tend to have even higher [Ba/Sr] than the r-II 
stars. The r+s stars are, however, the most overabundant in Ba 
relative to Sr of all groups, with a mean of [Ba/Sr] ~ 1.6. 

We note from Table|8]that there may also be a difference in 
the metallicity distribution of the stars within these groups. The 

r-II stars are centred on a metallicity of [Fe/H] 2.84 with a 

scatter on the order of 0.16 dex (see also Paper n). The r-I stars 
are found across the entire range of metallicities investigated in 
Papern, i.e. -3.4 < [Fe/H] < -1.5. The r+s stars from Table|8] 

are centred on [Fe/H] 2.55 with a scatter on the order of 

0.26 dex. Thus, the r+s stars seem on average to have a higher 
metallicity and a wider range of metallicity than the r-II stars. 

6. Discussion 

In this section we will discuss both the results for 
HE 0338-3945, and the r+s stars in general. We first sum- 
marise the abundances of HE 0338-3945 and then briefly dis- 
cuss the classification of the r+s stars, in the light of our survey 
of stars with n-capture enhancement in Sect. 15.41 Finally, we 
present a detailed discussion of the possible formation scenar- 
ios for the r+s stars, based on both our results for the abundance 
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Fig. 8. Plots showing the ratio of the different heavy n-capture 
i-elements (Ba, La, Ce) compared to the r-element Eu for stars 
of different classes. Symbols and data as in Fig. [7] 



pattern of HE 0338-3945 and our survey of r+s stars from the 
literature. 

6.1. The abundances of HE 0338-3945 

Abundances were derived for 33 elements and we esti- 
mated upper limits for an additional 6 elements for the star 



Fig. 9. Histograms showing the distribution of different classes 
of stars when comparing the heavy n-capture s-element Ba to 
the light n-capture ^-element Sr. Data as in Fig. [7] 



HE 0338-3945. We have confirmed the high content of the r- 
elements Eu, Gd, Dy, Er, and Tm with a mean of [r/Fe] ~ 2.10 
([Eu/Fe] ~ 1.94). The upper limits on Tb, Ho, Lu support an 
overall high r-element abundance, although the Ag upper limit 
is surprisingly low and not consistent with a scaled solar r- 
process distribution as normalised to Eu. This could, however, 
be an effect of over-ionisation of Ag I. We could detect neither 
Th nor U, as the spectral lines of these elements are weak and 
reside in heavily blended regions; only rather high upper limits 
were obtained. 

Although the star is also enhanced in ^-elements, there is 
a pronounced difference among the Is and the hs nuclei. The 
mean of the Is overabundance is [Sr, Y, Zr/Fe] = 0.92, and the 
hs [Ba, La, Ce/Fe] = 2.27. Obviously, the light ^-elements do 
not fit the solar s-element distribution of s-elements when nor- 
malised to Ba. However, they do fit the r-element distribution 
as normalised to Eu. This behaviour has been noted in other 
stars as well, such as HE 2148-1247 (Cohen et al. 120031 . The 
abundance of Pb is higher than expected from scaled solar s- 
element abundances as normalised to Ba, by about 1 dex. The 
elements that for the Sun were produced in significant amounts 
by both the r- and s-processes (e.g. Nd, Sm and Yb) are not 
unexpectedly overabundant compared to both the r- and the s- 
element solar distributions. 

The overabundance of C is confirmed, and we have derived 
an isotopic ratio of 12 C/ 13 C~ 10. Also N and O are overabun- 
dant in the star. Na to Ni seem normal for matter contami- 
nated by SNe of type n, the exception being Sc which is sig- 
nificantly enhanced compared to normal Population n stars (but 
see Sect. 14. 1> . The abundance of Cu is low relative to iron, and 
comparable with results for normal stars of similar metallicity 
(Mishenina et al. l2002l Bihain et al. l2004l 

The NLTE and 3D effects discussed in Sect. 14.11 may in 
many cases be severe. We lack estimates of corrections for 
most elements, and, moreover, the estimates we have are of- 
ten quite uncertain. It is thus risky to draw conclusions about 
nucleosynthesis sites by comparing abundances of different el- 
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Table 9. Definition of classes for n-capture-rich stars. In this 
article we define stars not falling under any of these classes as 
"normal". 



Class 


Constraints in [X/Fe] 


Other constraints 


r-1 


+0.3 < [Eu/Fe] < +1.0 


[Ba/Eu] < 0.0 


r-H 


[Eu/Fe] > +1.0 


[Ba/Eu] < 0.0 


r+s 


[Ba/Fe] > +1.0 
[Eu/Fe] > +1.0 


[Ba/Eu] > 0.0 


s 


[Ba/Fe] > +1.0 
[Eu/Fe] < +1.0 


[Ba/Eu] > 0.0 


Pb 


[Ba/Fe] > +1.0 


[Pb/Ba] > +1.0 



ements. One may hope that, e.g., the net effect of NLTE and 
3D could well be of the same order of magnitude for the Is- 
element Sr and the /w-element Ba, but the different character 
of the lines used may speak against such an optimistic view. 
Comparisons between different groups of elements, when car- 
ried out differentially relative to other stars with similar tem- 
peratures and gravities and overall metallicities may, however, 
be discussed with some confidence, in particular when we have 
some guidance from theoretical estimates from NLTE and 3D 
calculations. Unfortunately, 3D corrections for Eu are unavail- 
able at present, and these effects could perhaps be of some im- 
portance. 

6.2. The class r+s stars 

With the background of Sect. l5.4l we shall here discuss the clas- 
sification of r+s stars. We propose to use [Eu/Fe] as a measure 
for the r-elements, and [Ba/Fe] for the .s-elements. This choice 
of classification criteria, depending on single n-capture element 
abundances instead of e.g. all /i.s-elements, is motivated by sim- 
plicity and practicality. It could be of value to base a classifi- 
cation system on many criteria to suppress effects of errors and 
individual abundance fluctuations, but in practice both Eu and 
Ba have strong lines which are easy to observe in metal-poor 
stars and moreover, ambiguities are minimised. Also, this sim- 
ple classification could be transformed to a purely empirical 
system, independent of uncertainties in the abundance analy- 
sis. 

The criteria chosen for the r+s stars are presented in 
Table[9] It is natural to set the limit for .s-elements as [Ba/Fe] > 
+ 1.0, see Fig.0 but the limit for the r-elements may be de- 
bated. For simplicity we have made a choice analogous to that 
for the .s-elements, and require [Eu/Fe] > +1.0, although there 
is no natural separation between the r+s stars and the s stars 
in this respect. However, the distribution of the stars in the 
[Ba/Eu] versus [Ce/Eu] and [La/Eu] versus [Ce/Eu] diagrams 
(see Hig-inj suggests that there might be a natural separation of 
these classes of stars. As seen from the plots of Fig.[8]it is also 
natural to set the criteria [Ba/Eu] > 0.0 for the r+s stars. Note 
that we do not require the r+s stars to be rich in C, although all 
such stars known to us have [C/Fe] > +1.7 



The definition of the class r+s only partially overlaps with 
the class "r/s stars" defined by Beers & Christlieb (2005 1. The 
r/s criteria is set to 0.0 < [Ba/Eu] < +0.5, and thus also in- 
cludes stars enhanced in neither r- nor .s-elements. The r+s 
class, on the other hand, includes stars with [Ba/Eu] > +0.5. 

As seen in Table|9j we use the same criteria for r-I and r-II 
stars as Beers & Christlieb ( 2005 ), but this is not the case for s 
stars. Stars highly enhanced in .s-elements, only moderately in 
r-elements, and with [Ba/Eu] > 0.0 are called s stars here. This 
definition differs from the classification of Beers & Christlieb, 
who define s stars as having [Ba/Fe] > +1.0 and [Ba/Eu] > 
+0.5. 

Lead stars are denned as [Pb/hs] = [Pb/Ba, La, Ce] > +1.0 
by Van Eck et al. (2003 1. According to the same principles as 
discussed above, we propose to use only Ba as the measure for 
the /i.s-elements, and thus advocate the classification given in 
Table[9] We do not require the Pb stars to also be r+s stars, 
only that they are enhanced in Ba and Pb, and much more in 
the latter element. 

6.3. Scenarios for the formation of r+s stars 

The abundance patterns of stars enriched in both r- and s- 
elements may give important clues to the sites of nucleosyn- 
thesis for these elements. Many different scenarios have been 
proposed, most of them invoking two more or less indepen- 
dent processes, e.g. with one site for the production of the 
r-elements and one for the .s-elements. A basic question is 
whether such scenarios may be retained, in view of the ho- 
mogeneity of the r+s group, and its separation from the other 
classes of Population n stars. Another factor of significance is 
the frequency of r+s stars, relative to stars with more normal 
Population n composition. From Paper n we estimate that on 
the order of 1 % of the Population n stars are r+s stars, though 
we note this may be significantly underestimated due to the bias 
against CH strong stars. Here, these different r+s scenarios will 
be discussed. 

I. Radiative levitation: Peculiar abundance patterns may arise 
due to radiative levitation of elements with many absorption 
lines, while heavy elements with few lines may sink in the at- 
mosphere. This is seen in some Ap stars, which display enor- 
mous enhancements of the rare-earth elements. 

Cohen et al. (2003 1 discussed this possibility as an explana- 
tion for the r+s stars and dismissed it. The r+s stars do not show 
the same abundance patterns as the Ap stars. They presumably 
are not hot enough to produce any significant radiative levita- 
tion in their atmospheres. In particular, they have much deeper 
convective zones, efficiently mixing and diluting the surface 
material with the stellar interior. Note also that there are sev- 
eral r+s stars on the giant branch, where the convective zone 
contains most of the stellar mass. 

II. r-rich ISM and self-pollution: The star may have been 
formed out of r-enriched material from e.g. an early supernova, 
and later self-enriched its surface with .s-elements during He 
shell flashes in the AGB phase (as discussed by Hill et al. 2000 
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Cohen et al. 2003). We note that low-mass metal-poor stars 
may also be affected by dredge-up at the He flash (Fujimoto 
et al. l2TjtjuV However, both HE 0338-3945 and HE 2148-1247 
are located in the turn-off region and presumably have not 
passed the red-giant phase yet, so in the case of at least these 
r+s stars it seems possible to dismiss the hypothesis. 

III. Binary system out of r-rich ISM and AGB-pollution: The 

star could have been a secondary in a binary system, formed 
out of r-enriched material, and later polluted with C and s- 
elements by mass transfer from its AGB companion (discussed 
by e.g. Hill et al. l2000l Cohen et al. l2003l Ivans et al. 12005V 
We note in passing that metal-poor stars may well produce s- 
elements even without any Fe seed nuclei, according to Siess 
et al. d2002l 

The mass transfer from a companion AGB star is the sce- 
nario generally adopted for the formation of CH stars, known to 
be single-line spectroscopic binaries (McClure 1983 1984|and 
McClure & Woodsworth 1990 1. It is not clear, however, if all 
r+s stars also are CH stars, although all 17 r+s stars in Tabled 
have considerably enhanced carbon abundances, as noted in 
Sect. 16.21 Many CH stars are defined as s stars; the s stars are 
often quite C rich (see Table|8|l, while the r stars are generally 
less C rich. A relatively high fraction of the CH stars seem, 
however, to be r+s stars (e.g. Aoki et al. 2002). 

There are severe constraints set on the period and orbit for 
this mechanism to result in a C and s-enriched matter trans- 
fer, as the mass transfer should occur in the AGB phase where 
the ^-elements are produced and not (only) in the RGB phase 
(Jorissen & Boffin 1992). The transfer across the binary has 
to be efficient to reach the levels of contamination seen in r+s 
stars. Also the right amounts of elements should be transferred. 
A glance in Table[8] shows that the carbon enrichment is re- 
markably constant among the r+s stars. The low 12 C/ 13 C ratio 
of HE 0338-3945 suggests that a rather large fraction of the 
carbon (which presumably was produced by He burning in the 
primary) has been CNO-processed; N and O abundances also 
clearly suggest heavy processing but the dominant C relative to 
N and O puts constraints on this. There are obviously several 
steps that have to work, from the initially r-enriched material 
from which the star happened to form, via the formation of 
carbon, the formation of the ^-elements in the AGB inter-shell 
nucleosynthesis, the CNO processing of C and the pollution of 
the surface of the companion with sufficient amounts of matter. 

One may discuss whether this is a realistic scenario, in 
view of the many r+s stars discovered. In general, if a mech- 
anism with two independent steps is in action, one producing 
the r-enrichment and another the s-enrichment, and the two are 
stochastically independent, one would expect the probability to 
find a star affected by both mechanisms, p(r, s) to be the prod- 
uct of the probabilities of finding a star subjected to only one 
of the two, p(r) and p(s), respectively, i.e. p(r, s) = p(r) ■ p(s). 
Presently, in the HERES survey there are 8 r-H stars, two s 
stars and three r+s stars, among the 253 stars without strong 
CH lines analysed in Paper n. The CH-strong rejection criteria 
may lead to a strong bias, suppressing the number of s stars, and 
even r+s stars. In order for a two-step mechanism to produce 3 



r+s stars, however, a considerable fraction of all Population n 
stars should show s-star characteristics. In fact, we find from 
elementary statistical arguments that the probability of finding 
3 or more r+s stars out of 250 would be only 0.9% if the frac- 
tion of s stars would be 0.1 of all Population n stars. Even if 
the total fraction of s stars were 0.3, the chance to find 3 r+s 
stars would be only 22%. One should also note that in Paper n 
the r+s stars can only be detected when they are warm enough 
such that the CH lines are not so strong as to render the method 
used uncertain. Only a small fraction of the r-II stars found in 
the HERES survey are so warm, which suggests that the ratio 
of stars with r+s abundances relative to r-II stars is consider- 
ably higher than suggested by the HERES results. Thus, the 
probability that an r+s star is just a star which happened to be- 
come r-enriched, e.g. due to initial inhomogeneities in the ISM 
or a later nearby SN, and independently of that had an AGB 
enrichment from a companion, is negligible and should be dis- 
missed. Still another argument against such a two-step scenario 
is the chemical homogeneity within the r+s group and with a 
clear separation from r-II stars and s stars. The scenario is also 
problematic as it is not commonly accepted that the ISM was 
sufficiently inhomogeneous (e.g. Qian & Wasserburg 2001) to 
contain such great s- and r-element overabundances relative to 
Fe at [Fe/H]> -3. 

A possible modification of this scenario would be that 
the formation of the binary system is triggered by a su- 
pernova which also provides the r-elements (discussed by 
Gallino e t al. 120051 and Ivans et al. [2005 1. Vanhala & 
Cameron d!998i have performed simulations of triggered star 
formation due to shock waves impacting molecular cloud 
cores, suggesting that weakly evolved cores may fragment dur- 
ing collapse and form low-mass binaries. Such a scenario is ap- 
pealing since the r- and ^-enrichments are no longer stochas- 
tically independent and thus the high frequency of r+s stars 
might be explained. However, it is presently unclear whether 
this scenario could be so prevalent as to explain the observed 
frequency. 

Cohen et al. (20031 argue for a ^-enhancement due to 
an AGB contamination, although they prefer scenario V be- 
low. They point out that several of the stars we now de- 
fine as r+s stars are binaries (CS 29526-110, HE 2148-1247, 
HE 0024-2523, LP 625-44), although we cannot confirm this 
for HE 0338-3945, nor for several other r+s stars. Preston & 
Sneden (2001 ) report on r+s stars not exhibiting radial veloc- 
ity variations exceeding 0.5 km/s over an 8 year period. This 
discussion may also be applied to scenarios IV, V and VI. 

IV. Triple system with SN- and AGB-pollution: The star could 
have been a tertiary in a triple system, in which the primary 
exploded as a r-element producing supernova, and the other 
companion evolved into an AGB star dumping s-rich material 
on the least massive star. 

This scenario is also discussed and dismissed by Cohen 
et al. (2003 1 as being not very plausible. The likelihood that 
a star would survive a close SN explosion, that the secondary 
would not drift away, and that the secondary could subse- 
quently transfer processed matter onto the star at the right time 
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does not seem very great. If such events do happen, it is not 
probable that they would be common enough or sufficiently 
constrained to create so many r+s stars with similar abundance 
patterns, notably with similar relative amounts of C and of r- 
and s-elements. 

V. Binary system with AGB- and 1.5 SN-pollution: The star 
could be a secondary in a binary system where the companion 
first contributed s-elements as an AGB star, and later exploded 
as a "Type 1.5" supernova (Zijlstra 2004 Wanajo et al. 2005 1, 
producing the r-elements. The SNe of type 1.5 are proposed to 
be more common among metal-poor stars than for Population i, 
due to the strong metallicity dependence of the mass loss dur- 
ing the AGB phase leading to a different initial-final-mass rela- 
tion for low-metallicity stars. Thus, intermediate-mass metal- 
poor stars may end up with higher final masses after the AGB 
phase, and more easily reach the Chandrasekhar mass and ex- 
plode as a Type 1.5 supernova. 

VI. Binary system with AGB- and AlC-pollution: The star 
could be a secondary in a binary system, in which the com- 
panion first contributed s-elements as an AGB star and then 
evolved to a white dwarf. Later, the primary might have un- 
dergone an accretion-induced collapse (AIC) in the white 
dwarf stage, triggered by mass transfer in the other direc- 
tion across the system. The collapse created a neutrino wind 
and r-elements which were transferred across the binary 
and contaminated the surface of the secondary star (Qian & 
Wasserburg 2003 Cohen et al. 2003 1. Cohen et al. also noted 
that the system might disrupt during the AIC event, which 
would explain an apparently single r+s star. 

This hypothesis is physically uncertain due to uncertainties 
in neutrino and neutron-star physics (Qian & Woosley 1996 1, 
and it is thus uncertain whether it works at all. It is also unclear 
how the white dwarf would accrete matter from the secondary 
star, which in several cases has been found not to be evolved 
beyond the turn-off phase (HE 0338-3945, HE 2148-1247). 
One may question how probable is it that the white dwarf is 
close enough for mass transfer, in view of the mass-loss it has 
experienced which increases the distance between the compo- 
nents. 

VII. Binary system with only AGB-pollution: The star might 
be a secondary in a metal-poor binary system with an AGB 
companion producing n-capture elements in a hypothetical 
high neutron density s-process. The relatively high flux of 
neutrons has been proposed to not only produce "normal" s- 
elements such as Ba and La, but Eu which otherwise is mainly 
assumed to be produced in the r-process with much higher neu- 
tron densities. 

This hypothesis has been discussed by Aoki et al. (|2002 1, 
Cohen et al. (2003 ) and Johnson & Bolte JIUOH . Aoki et al. 
and Cohen et al. found some support from the calculations 
by Goriely & Mowlavi (2000 1. Cohen et al. found a ratio of 
Ba/Eu ~ 100 for the star HE 2148-1247 (we similarly find 
Ba/Eu ~ 115 for HE 0338-3945), and Goriely & Mowlavi 



predicted a surface abundance ratio of 105 for their most metal- 
poor (Z=0.001) AGB-star model with dredge up from a long 
sequence of He shell flashes. This was, however, the resulting 
ratio from a relatively mild enhancement (0.48 dex in Ba, 0.07 
in Eu) of material with initially solar Ba/Fe and Eu/Fe ratios. 
The dredge-up material in the model of Goriely & Mowlavi has 
a Ba/Eu ratio of about 500. Anyhow, Cohen et al. dismissed the 
scenario since they found the predicted absolute abundances of 
the produced elements to be quite insufficient to explain the ob- 
servations. Johnson & Bolte, on the other hand, found evidence 
in the r-element ratios (Eu/Tb, Eu/Dy, Th/Eu) which seem dif- 
ficult to explain with any combination of the normal r- and 
s-processes. They suggested that another form of s-process is 
at work, and found some support for this in parametrised calcu- 
lations made by Malaney ( 1987 1. In Malaney's Table II, yields 
of the s-process acting at the unusually high neutron density 
of 10 12 cirT 3 are presented, which reduced the Ba/Eu ratios to 
150 - 200. Adopting a neutron exposure of tq = 0.05 Johnson 
& Bolte found the predictions of Malaney to agree reasonably 
well with the observed La/Eu ratio for the r+s star CS 31062- 
050, and a generally good overall fit of the abundance pattern of 
that star from La to Hf. They noted, however, that the observa- 
tions of the 151 Eu/ 153 Eu ratio of this star by Aoki et al. (|2002i 
seem to indicate that the neutron density was rather around the 
more normal value of 10 8 . 

In order to study the consequences of a high neutron- 
density s process we set up a small program to calculate s- 
process yields in the interval from La to Gd, simultaneously 
solving the rate equations for 88 nuclides in this interval. The 
n-capture cross sections were adopted from Bao et al. (|2000l 
with some modifications such as the new data from Best 
et al. (2001 1. The /3-decay rates and the electron-capture rates 
were from Takahashi & Yokoi (1987). In stationary-flow cal- 
culations we found the abundances of Eu relative to La, Ce, Pr, 
Nd, Sm and Gd to stay significantly below the observed values 
and vary rather little when the neutron density increased from 
10 8 to 10 14 cirT 3 , with some tendencies to peak around 10 12 ; 
only Eu/Sm and Eu/Pr came then close to the observed val- 
ues. For the time-dependent case we varied neutron exposures 
and temperatures, but again in this case we could not recover 
the high observed Eu/X ratios for X= La, Nd, Ce, or Gd with 
any combination of parameters tried. Although these results are 
preliminary, and a number of nuclear rates are still rather un- 
certain, we find it less probable that a high-neutron density s- 
process would be responsible for the characteristics of the r+s 
stars. 

VIII. A relation to blue stragglers? The fraction of r+s stars 
in Table [8] that have T e g > 6000 K is fairly high, as com- 
pared with the fraction of hot s stars and r-II stars. This may re- 
flect small-number statistics or selection effects - it is certainly 
more difficult to detect an r-II star than a r+s star at the turn-off 
point (TOP) due to the fact that the former, at least in our sam- 
ple, tend to have both lower metallicities and lower [Eu/Fe]. 
Furthermore, since r+s stars typically have [C/Fe] > 2.0 (see 
Table|8|l, the cooler r+s might have been rejected from the sam- 
ple of Paper ii due to their strong CH lines. The dominance of 
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detections of r+s stars at the TOP might also point at some 
relation between the r+s stars and the blue stragglers in glob- 
ular clusters and in the field. These objects, lying above and 
blueward of the turn-off in the colour-magnitude diagram of 
clusters, are generally believed to be the result of close-binary 
evolution with mass transfer, of stellar collisions or of mix- 
ing of the stellar interior (for a review, see, e.g. Trimble 1993 1. 
Sneden et al. (2003b I recently analysed spectra of six field blue 
stragglers with [Fe/H] < -2. Three of them were spectroscopic 
binaries and were found to have enhanced carbon and Sr and 
Ba. This was interpreted as the result of mass transfer from 
companions in the AGB stage. One of these stars was remark- 
ably rich in Pb and had in fact a high, though uncertain, Eu 
abundance - that star is CS 29497-030 and is included as an 
r+s star in our Table [8] The idea that sub-giant CH stars may 
be related to blue stragglers was already put forward by Luck 
& Bond ( 1991 1. However, in order for this scenario to apply to 
the r+s stars, and not only as a provider of the ^-elements, an 
idea which is already included in the scenarios III- VI above, 
one has to invoke a production mechanism for the r-elements. 
One might hope that the alternative origin of blue stragglers 
as the result of stellar collisions might provide not only over- 
all mixing of CNO and He (e.g. Sills et al. 12005b but possi- 
bly also heavy neutron fluxes at head-on collisions, leading to 
formation of Eu and other r-elements. It is, however, highly 
uncertain whether such a speculative mechanism, presumably 
acting in the centre of dense globular clusters, could contribute 
enough r+s stars in the field. 

IX. Origin from a common cloud of r+s stars: In view of the 
similarity in chemical composition of the r+s stars it is not 
inconceivable that they were formed out of gas in a common 
cloud, and even were accreted by the Galaxy at a later stage. 
To explore this hypothesis we have calculated space velocities 
for the stars. In the absence of accurate distances for but one 
star (CS 31062-012), we used photometric distances estimated 
by reading off the colour-magnitude diagram of the metal-poor 
globular cluster M92 (Ruelas-Mayorga & Sanchez 2005 1 at the 
appropriate reddening-corrected B - V. Alternatively, we have 
compared to theoretical isochrones of Kim et al. (|2002 ) with 
[Fe/H] = -2.5 and an age of 12 Gyr in the r e ff-logg diagram 
and derived the absolute magnitudes from that. For the four 
stars with well known kinematic data we find space velocities 
with a spread of approximately 100 km/s in the U, V and W 
components, and conclude from this that there is no suggestion 
of a common space motion. 

7. Conclusions 

In this paper we have analysed a high-quality spectrum of the 
r- and ^-element rich star HE 0338-3945, deriving abundances 
for 33 elements and upper limits for an additional 6 elements. 
We found the abundances of this star to be very similar to 
those of HE 2148-1247, which is supported by an homoge- 
neous analysis of the two spectra. In fact the abundance pat- 
terns among known r+s stars in the literature seem to be quite 
similar considering the differences between analyses. Based on 



our results for HE 0338-3945 and our survey of other r+s stars, 
we have discussed a range of scenarios for the formation of the 
r+s stars. Some scenarios have considerable merit, but all have 
also drawbacks or problems. Scenario I is already dismissed. 
Scenarios II and III need the ISM to be highly inhomogeneous 
to produce such high r-process abundances as detected in these 
stars. We note that the [Eu/Fe] ratios for the r+s stars are sys- 
tematically higher than those of the r-II stars which might sug- 
gest that a process with two independent steps, one of which is 
also responsible for the r-II stars is not very probable. However, 
this is uncertain due to the possibly significant contribution of 
Eu by the .s-process, and the possible relative effects of con- 
vective mixing between the r+s stars, which are predominantly 
TOP stars, and the r-II stars, which are all giants (see also dis- 
cussion below). Scenario II is also improbable as it requires 
all r+s stars to have passed at least the RGB phase. Scenario 
III presently seems unlikely in view of the frequency of r+s 
stars and their apparent chemical homogeneity. The discussed 
modification of scenario III to include supernova-triggered bi- 
nary formation seems plausible, though it is presently unclear 
whether such events can provide the observed frequency of 
r+s stars. Scenario IV is improbable as it requires strict stel- 
lar orbit constraints to produce the chemical homogeneity of 
the r+s class. Scenarios V and VI may be considered, but as 
they depend heavily on physical processes and parameters that 
are poorly understood at present, these scenarios at best have to 
be regarded as uncertain. The observed chemical uniformity is 
not an obvious consequence of them. Scenario VII still needs 
a site which can produce the high neutron densities at suitable 
exposures; the scenario also has to be supported by renewed 
and more realistic ^-process calculations. It is also not clear 
how well this scenario would explain the chemical homogene- 
ity. A probable merit of it is, however, that the observed great 
but varying amounts of Pb may result naturally from the high 
neutron fluxes. The absence of observed radial velocity varia- 
tions for several stars (see Preston & Sneden 120011 and Aoki 
et al. 2002 1 is a problem, although not very severe as yet, for 
all scenarios III- VIII. Scenario VIII seems hypothetical and it is 
also questionable whether it is efficient enough. The available 
astrometric data, though admittedly meager, does not support 
Scenario IX. 

The majority of the r+s stars in Tableware TOP stars, while 
all the r-II stars are giants. One should note the great difference 
between the TOP stars and the red giants as regards two im- 
portant aspects: First, the convective zone of the TOP stars is 
less than 5% of the stellar mass, while it is typically 75% of 
the mass for the giants. Thus, the dilution of enriched mate- 
rial transferred to the stellar surface of a TOP star is at least a 
factor of 10 less than for a giant. So, if the enhancements of 
n-capture elements are due to mass transfer from a companion, 
the amount transferred has to be at least one order of magnitude 
greater for the r-II stars than for most of the TOP stars. Second, 
the space volume surveyed in apparent-magnitude limited sur- 
veys is at least a factor of 200 greater for the giants than for 
the r+s stars. This suggests that the space density of r+s stars 
is considerably higher than that of r-II stars. In spite of their 
scarcity, it is interesting that no r-II stars have been found at 
the TOP, since a much smaller transfer of r-elements would 
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raise the atmospheric abundances to considerable values. More 
statistically controlled surveys are desirable. 

Future detailed studies, with accurate analysis of all known 
and suspected stars, are worthwhile. Spectra with high S /N and 
high resolution, also in the the UV, are important to obtain bet- 
ter limits on Ag, Th and U and to get results for r-elements that 
have lines in the UV such as Os and Ir. Attempts to measure Eu 
isotopic ratios for these stars should be continued. The analy- 
sis of certain important elements, such as Eu and Ag, should 
be scrutinised for NLTE and 3D effects. More detailed calcula- 
tions of the i-process at high neutron densities should be car- 
ried out. Finally, all known r+s stars should be monitored in an 
attempt to disclose their possible binarity. 
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Appendix A: Line List 

The complete line list is presented in Table 2, which is avail- 
able only electronically. Here we briefly comment on line se- 
lection issues and data sources for each element, particularly 
for the oscillator strengths, hfs and isotopic splitting. Other 
data such as wavelengths, excitation potentials were collected 
from various sources, often the Vienna Atomic Line Database 
(VALD, Kupka et al. ll999l . the NIST compilations or Sneden 
etal. (2003a}. 

Carbon: The abundance of carbon was derived from the A- 
X bands at 43 10^3 13 A (0,0) and 4362^1367 A( 1,1), along 
with two B 2 S-X 2 n (1,0) CH lines at 3638 A and 3661 A. The 
band regions used in this programme are the same as those used 
by Paper ii. The two single lines were selected from the list 
of J0rgensen et al. 1 1996) as they are unblended, of moderate 
strength, and have reasonable S /N in the observed spectrum. 
The C abundance was derived from 12 CH features assuming an 
isotopic ratio of 12 C/ 13 C = °o as in Paper n. 

Nitrogen: The nitrogen abundance was derived using the 
B2-X2 bands of CN at 3872 A (1,1) and 3882 A (0,0). The line 
data was extracted from a list made available by Plez (1998, 
private communication) and is described by Hill et al. (2002 1. 
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Oxygen: The data for the oxygen lines were compiled from 
VALD, except for the logg/ values which were astrophysical 
values taken from Jonsell et al. (2005). The differences from 
the NIST log gf values in VALD were less than 0.02 dex. 

Sodium: For two of the three Na i lines the /-values of Wiese 
& Martin J1980I and hfs from McWilliam et al. dl995> were 
used. The /-value for the 8183 A line was taken from VALD. 
No hfs data were available for this line, and as it is weak this 
should not be important. 

Magnesium: A total of 7 lines of Mgi were used. Data for 
the lines at 3829, 4571, and 4703 A were taken from Paper n. 
Due to contamination by CH, Balmer lines, and some resid- 
ual reduction artefacts not all lines from Paper n could be used. 
Oscillator strengths for the lines at 5173, 5184, and 5528 A 
originate from Wiese & Martin ( 1980 1. For the 5711 A line 
the astrophysical logg/-value from Jonsell et al. d2005i was 
adopted. 

Aluminium: The Al i resonance line at 3961 A was used for the 
analysis. The hfs was taken from McWilliam et al. (1995 1. 

Calcium: The /-values for the 10 Cai lines were taken from 
Wiese & Martin (fT980l . Smith (TrWA . Jonsell et al. J20051 and 
VALD. 

Scandium: The abundance of scandium was determined from 
5 Sen lines, using /-values from Lawler & Dakin (1989). Hfs 
for the lines 4247 and 4415 A was taken from McWilliam 
et al. (I1995> . The remaining lines are very weak and hfs was 
not considered. 

Titanium: A total of 11 lines of Tii and 13 of Tin were 
analysed. The oscillator strengths were gathered from various 
sources: Blackwell et al. dl982a> . Blackwell et al. d!982b> . 
Fuhr & Wiese JT996I ). Grevesse et al. JT9891 . Martin 
et al. ( fT988l . Pickering et al. d200Tl . Ryabchikova et al. (fT994l ) 
and VALD. 

Vanadium: The vanadium abundance is based on one Vi 
and 4 Vn lines. The /-value data were taken from Fuhr & 
Wiese (tT996V 

Chromium: The /-values for the 6 Cri and 2 Cm lines were 
taken from Fuhr & Wiese ( 1996 ) and Martin et al. ©88 ). 

Manganese: For the two lines of Mn i we adopted the log gf 
data from Booth et al. ( 1984 1, and hfs was computed using the 
data compiled in Lefebvre et al. ( 2003 1. For the three Mn n lines 
we used data from Kling & Griesmann ( 2000 ), with hfs from 
Holt et al. JT9991 



Iron: A total of 61 Fei lines and 8 Fen lines were anal- 
ysed. The neutral line oscillator strengths are taken from 
O' Brian et al. fTWH . Fuhr& Wiese ( fT996l . Bard et al. JT99TI . 
Bridges & Kornblith dl974i . and VALD. For the Fen lines 
/-values were adopted from Biemont et al. ( 119911 ). Fuhr & 
Wiese (fT996l . Heise & Kock ©90), Kroll & Kock tTW7\ . 
Moity fT983l Schnabel et al. i2004t . and VALD. 

Cobalt: The cobalt abundance was derived from 4 Coi lines 
using /-values from Cardon et al. ( 1982 1 and Nitz et al. ©99 ), 
and hfs data from Pickering ( 1996). 

Nickel: The 8 lines of Nil have oscillator strengths from 
Black well e t al. (1989}, Fuhr et al . (IT9881 . Fuhr & 
Wiese JT9961 . Huber & Sandeman <fl980t . and VALD. 

Copper: We analysed two lines of Cui, both in the UV. The 
log gf data were taken from Fuhr & Wiese ( 1996 ), and hfs data 
from Kurucz ( 1995 ). 

Strontium: One Sri line and 4 Srn lines were analysed. 
Hyperfine and isotopic splitting are accounted for in the 4077 
and 4215 A lines (see discussion in Papern). We did not find 
hyperfine structure data for the remaining lines; however, as 
they are weak, this should not be important. The /-values 
were gathered from Migdalek & Baylis (1987), Pinnington 
et al. (IT9951 and Sneden et al. d2003a> . 

Yttrium: The /-values for the 11 lines of Yn are from 
Hannaford et al. (tT982l 

Zirconium: All but two of the 16 lines of Zrn were analysed 
with /-values from Biemont et al. d!981> . The remaining lines 
at 3356 and 4019 A have /-values from VALD. 

Silver: The upper limit to the silver abundance was derived 
from two Ag i lines. The oscillator strengths were taken from 
from Fuhr & Wiese d!996i . The hfs data was taken from Ross 
& Aller JT9721 . 

Barium: The barium abundance is derived using three Ban 
lines. The log gf data came from Gallagher ( 1967 1 and VALD. 
The isotopic splitting and hfs were taken from Villemoes 
et al. ( 1993 ) for the two weak lines, noting that data for only the 
lower levels are available. For the 4554 A line, hfs and isotopic 
splitting from McWilliam ( 1998 ) was adopted. In all cases, the 
solar isotopic composition from Anders & Grevesse ( 1989 ) was 
assumed. We note that if a pure r-process isotopic composition 
is assumed we find an abundance only 0.03 dex lower. 

Lanthanum: The 8 lines of Lan were analysed employing os- 
cillator strengths and hfs data from Lawler et al. (12001a). For 
the lines at 4921 and 4922 A only the hfs of the lower levels 
was considered, as no data were available for the upper levels. 
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Cerium: Abundances were derived from 15 lines of Cen. 
The /-values were gathered from various sources, Palmeri 
et al. (HHP, Sneden et al. {TH), Gratton & Sneden ( 173331 . 
and VALD. 

Praseodymium: The logg/ values for the 4 lines of Prn were 
taken from Ivarsson et al. ( 1200 1 i . 

Neodymium: A total of 34 Ndn lines were analysed. The os- 
cillator strengths were taken from Den Hartog et al. (2003 1. 

Samarium: The 4 lines of Sm n were analysed using /-values 
from Biemont et al. (I1989> . 

Europium: The abundance of europium is based on 3 Eun 
lines, with /-values, hfs and isotopic splitting taken from 
Lawler et al. (2001ci. Solar r-process isotopic fractions were 
assumed. 



H af n i u m : We have one line of Hf n with the log gf taken from 
VALD. 

Lead: The abundance is derived from two Pbi lines. The os- 
cillator strengths were taken from Biemont et al. ([2000 » and 
the hfs and isotopic splitting data from Manning et al. |l950l. 
Solar isotopic ratios from Anders & Grevesse ( 1989 1 have been 
assumed. 

Thorium: An upper limit for the thorium abundance is de- 
rived from the Th n line at 4019. 129 A employing the oscillator 
strength from Nilsson et al. ( 2002b I. 

Uranium: The /-value for the uranium line is taken from 
Nilsson et al. H2002a> . 



Gadolinium: The /-values for the 8 Gdn lines were collected 
from Sneden et al. ( 2003a i; the paper of Cowan et al. (|2002 1 
should be consulted for details. 



Terbium: The /-value for the single line of Tbn was taken 
from Lawler et al. |2001di. 

Dysprosium: The /-values for 15 of the 16 Dyn lines anal- 
ysed were taken from Wickliffe et al. ( I2000> . For the line at 
3944 A the /-value is from from Biemont & Lowe (1993 1. 

Holmium: The single Hon line at 4152 A has atomic data in- 
cluding hfs taken from Sneden et al. (2 003a> . 

Erbium: The /-values for the 3 Ern lines were taken from 
Musiol & Labuz ( fT983l and Kurucz (1995}. 

Thulium: The 4 Tmn lines have oscillator strengths from 
Kurucz (1995 1. According to Sneden et al. (1996 1 these are 
rescaled laboratory data from Corliss & Bozman I1962> . 

Ytterbium: The /-values used for our two Yb n lines are from 
Sneden et al. ( 2003a I. They have renormalised theoretical os- 
cillator strengths from Biemont et al. (1998 1 using lifetimes 
from Pinnington et al. ( 1997 1. The isotopic and hyperfine split- 
ting data are taken from Martensson-Pendril et al. ( 1994 1, and 
the solar isotopic ratios from Anders & Grevesse ( 1989 1 have 
been assumed. 



Lutetium: For the single Lun line we adopted the log gf and 
hfs data from Sneden et al. ( 2003a I. The oscillator strength is 
from Quinet et al. (1999 1, renormalised using lifetimes from 
Fedchaket al. (2002). 



